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Abstract 

This paper presents a new approach on passive identification of elastic propagation media. 
As passive identification relies upon noise field correlation, an original and perhaps more natural 
approach is developed which consists in considering the Green correlation. It is shown that 
Green correlation contains all medium parameters and provides an appealing alternative to 
the classical Green function estimation through a Ward identity. Our formalism allows to 
extend classical scalar passive identification models to vectorial ones and to take into account 
realistic dissipation models. This approach is applied to acoustic and solid waves with viscous 
damping. 

1 Introduction 

Passive identification of a propagation medium consists in retrieving medium parameters, 
by using uncontrolled noise fluctuations [14]. Such an idea has long been pursued in acoustics 
([10], [14]) and seismology ([3], [15]) and gave rise to numerous applications and experimental 
validations ([13], [TTJ, [6]). 

Preceding studies ([10], [3], [13], [TTJ, [6]) rely upon the estimation of the Green function of 
the medium. Such estimation is made possible by exploiting a Ward identity [16], which relates 
the noise correlation function to the Green function ([10], [14], [5]). The fundamental role of 
dissipation in Ward identity was outlined in [6], were dissipation is assumed to be constant; 
however, a constant dissipation model is hardly acceptable from a physical point of view ([8], 
|12j). and needs to be further discussed. 

The contribution of the paper is twofold. Firstly, an alternative approach to passive iden- 
tification is proposed, which is based upon the Green correlation function. Green correlation 
is introduced as the correlation of a propagated white noise field [7]. We show that Green 
correlation contains all physical parameters that need to be identified for a complete character- 
ization of the propagation medium. The motivation for introducing Green correlation comes 
from classical system identification theory, where noise based identification relies strongly on 
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the transformation of second order statistics through linear systems [9] . 

Secondly, we highlight the role of Green correlation in the framework of elastic waves with 
viscous damping. Green function and Green correlation are explicitly formulated in the time 
and space Fourier formalism. Ward identities are also developed. Furthermore, for low viscous 
damping and far-field assumptions, a new formulation of the Ward identity involving third 
order time derivative, is derived. 

The organization of the paper is the following. In section 2, we introduce passive iden- 
tification through a linear system approach. We recall the definition of the Green function 
and its role in medium identification. Cross-correlation of random fields is recalled, and white 
noise notion is introduced. Then, we define the Green correlation, we show its role in passive 
identification and we discuss its relation with the Green function through Ward identities. 

In section 3, solid waves equation with viscous damping is presented. Dispersion matrix 
and relation of dispersion are introduced. Acoustic waves propagation is seen as a particular 
case of solid propagation study. 

Solid and acoustic Green functions are computed in section 4. We highlight their role 
in solid and acoustic medium identification, respectively. Approximation of far-field and low 
attenuation is also considered. This particular, but realistic, case allows to derive explicit ex- 
pressions. 

Solid and acoustic Green correlation are computed in section 5. The role of Green correla- 
tion in passive identification is emphasized in those practical cases. Ward identities are derived 
and compared to existing ones for a constant damping model. 

2 Green function and Green correlation of a linear 
propagation medium 

In this section, we adopt a linear system approach to describe medium identification. We 
introduce the Green function and we relate it to the classical impulsional response of linear 
systems. White noise is introduced and its existence is discussed. From white noise, we define 
the Green correlation and we highlight its natural role in passive identification. 

2.1 Medium and fields: a system approach. 

We denote by u(t, x) the value of field u at time t and position x. When u has one compo- 
nent it is said scalar, otherwise it is said vectorial. 

A propagation medium can be seen as a system where the source field f is the input and 
the generated field u is the output of the system. Those two fields are related by a relation of 
the following type: 

u=<9(f) (1) 

where operator & is a representation of medium properties: physical laws, boundary and initial 
conditions. The medium is said linear when is a linear operator. Only linear medium are 
considered in this paper. 

With that formalism, medium identification can be seen as system identification where 
parameters are physical, like attenuation or propagation speed, and geometrical, like distance 
or orientation between sensors. 
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2.2 Green function of a linear medium. 



A linear medium X satisfies the superposition theorem i.e. the value u(t,x) of the generated 
field can be seen as the superposition of all contributions of elementary sources f [if , x[)dt' dx[ 
emitted at time t' during dt' period in the volume centered in x' and of dimensions da:', for all 
times t' and points x'. Mathematically, this can be written u(t, x) — f R .. x G(t, x, t\ x')f(t', x')dt'dx' 
and simplified by introducing the generalised convolution £g>T,s as: 



G is called the Green function of the medium as it is the kernel of the operator appearing 
in equation ([TJ. All medium mechanical parameters are contained in its expression. This high- 
lights the importance of retrieving the Green function in medium identification. 

Physically, the i-th column of G corresponds to the medium response of a spatio-temporal 
impulsion directed by the i-th axis of the reference. According to this interpretation, G is 
sometimes called " impulsional response" of the medium in reference to the classical impulsional 
response of a linear system [9]. Then, active identification consists in emitting spatio-temporal 
impulsions to retrieve the Green function and then to estimate model parameters [9]. 

We consider only time-shift invariant media. This property implies that G depends only 
on the times difference appearing in its parameters i.e. we can do the following substitution: 
G{t,x,t',x') o G(t - t',x,x'). 

2.3 Cross-correlation of stochastic fields. 

In passive identification, source fields are not controlled. The principle relies on record- 
ing noise sources and using their statistical properties to retrieve medium parameters. With 
stochastic source fields, the analysis has to be performed from the cross-correlation of the 
generated field u defined as: 



where E and T are the expectation operator and the transposition operator, respectively. 

We consider only stationary fields, this assumption is not a strong constraint in practice. 
In that case, cross-correlation depends only on the times difference appearing in its parameters 
i.e. we can do the following substitution: C(t,x,t' ,x') <-> C(t — t',x,x'). Furthermore, when 
fields are also considered ergodic in time, cross-correlation can be computed using the following 
formula: 



This formula is fundamental to approximate the cross-correlation of a field recorded by an 
array of sensors. 

2.4 White noise. 

By definition, a white noise is a field which value of a given time, position and direction is 
uncorrelated to any other value taken at all other times, positions and directions. Mathemati- 
cally, the cross-correlation of a such field f is a spatio-temporal isotropic impulsion: 



u = G <g>T,S f 



(2) 




(3) 




(4) 



C { (t,x, t' , x) — 5(t,t')5(x,x')I 



(5) 
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where p is the number of components of f , I is the p x p identity matrix and 5 is the Dirac 
distribution. 

A white noise has no physical reality in a sense that it has an infinite power. However, 
in practice the temporal whiteness is only needed in a limited frequency band. This latter is 
defined by the used instrumentation. The classical approach [3] to justify that ambient noise 
converges to a source with a spatio-isotropic whiteness consists to see the medium as a chaotic 
dynamical system. Then, according to equipartition theorem, it exists a time after which a 
coherent source snared in the medium becomes spatially and isotropically white. This time, 
called mixing time, depends on the frequency band, medium geometry and heterogeneity. With 
those considerations, cross-correlation of ambient noise is stacked during a sufficient long time 
in order to obtain a contribution of an approximated spatio-isotropic noise [3J, [BJ. 

2.5 Green correlation. 

As the Green function is the field generated by a spatio-temporal isotropic source i.e. 
f (t, x) = 5(t, t')8(x, , we define by analogy the Green correlation C by the cross-correlation 
of a field generated by a white noise source. This function, introduced in [7], plays by definition 
a fundamental role in passive identification. 

We can precise Green correlation expression using definition ((3J and equation ([2]). Indeed, 
for every generated field u, C can be expressed as: 

C u =g® T ,sC ( ® T| sG" (6) 

where G~ (t,x,x') := G(— t, x, x') T . It is important to note that to establish j6]), we use the 
property: G(t,x,x') = G(t, x , x), true for all times t and all couples of positions (x,x'), 
according to spatial reciprocity. Equation (J6j) is the "order two" version of equation ((2}. 
When the source is a white noise, we obtain by using equation |[6j) an expression of the Green 
correlation: 

C:=G® T ,sGr (7) 

This shows the fundamental importance of the Green correlation in medium identification when 
statistical properties of ambient sources are taking into account. Equations ([6j and (J7J show 
that "perfect" white noise is for passive identification what "perfect" impulsion is for active 
identification. 

Green correlation does not appear in the literature as a fundamental field to retrieve in 
order to estimate medium parameters. Generally, the cross-correlation of a field generated by 
a white noise through a propagation medium is directly related to the Green function by the 
authors. This relation is called Ward identity and relates actually the Green correlation to the 
Green function. 

3 Elastic propagation with viscous damping. 

In this section, we recall acoustic and solid propagation with viscous damping equations 
|12| . Classical results are derived with a vectorial formalism which is useful to describe solid 
waves propagation, in particular the coupling between the P-waves and the S- waves. Dispersion 
matrix and relation of dispersion are presented in order to compute the acoustic and solid Green 
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function, the acoustic and solid Green correlation and Ward identities. Low attenuation case 
is also studied because it provides physical interpretations. 



3.1 Solid waves equation with viscous damping. 

We consider an elastic, homogeneous, isotropic and linear solid medium. Let (A, fi) be the 
Lame parameters expressed in N.m~ 2 , (x,v) be the viscous damping parameters for P- waves 
and S- waves, respectively, expressed in N.s.m~ 2 , and p be the density of the medium expressed 
in kg.m" 3 . Let f be a 3-components causal spatio-temporal displacement source and u the 
3-components displacement field. From Newton theorem and Hooke law, we obtain the solid 
waves equation with viscous damping [12] : 



I3 is the identity operator of 3-components fields, 



u = f (8) 



C := — - — VV T — — AI3 (Propagation operator) 
P P 

V := - X + r, X7X7 T - -AJ 3 (Dissipation operator) 
P P 

where A is the Laplacian operator and V is the gradient operator. 



According to the Helmholtz- Hodge theorem [5] , each displacement field u satisfying equation 
{HJ when there is no source i.e. f = 0, can be decomposed asu = V(/>p + V A ^> , respectively, 
where A is the cross product. <f>p is a scalar field, and, ip is a two-components field satisfying: 

4>p = (9) 

t s =° ( 10 ) 

with vp := ^j^, v% := E , a.% := 2t ^ 2 and a| := 2 . This shows that a solid wave is the 
contribution of two modes satisfying the classical acoustic waves equation with viscous damping 
|12j . <j)p is a pressure wave, called P-wave, of speed vp and attenuation ap. ijj is a shear 
waves, called 5- wave, of speed vs and attenuation as- 



9 2 9 A 2 A 

__ Qp _ A _ VpA 



° T 2 A 2 A 



3.2 Dispersion matrix and relation of dispersion. 

For a deterministic field u, we denoted by: 

u(w,fc) := / n(^)-£) e_1 ^ t_ ~ ~^dtdx (11) 

Jmxm 3 

its Fourier transform in the (w, fc)-domain where ui and k are the frequency variables associated 
with t and x, respectively. 

Applying time and space Fourier transform to equation (|8]), we obtain the algebraic relation: 

S(w,fc)u(w,fe) = ffw.fc) (12) 
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where is the dispersion matrix: 

:= (v%(u) - v|(cj))fcfc T + (v|(o;)A: 2 - w 2 )I| (13) 

with: Vp(oj) := up + iujap, v|(w) := u| +iwa|, k 2 := fc T fc and I3 is the 3x3 identity matrix. 

The dispersion matrix is fundamental as it plays the same role as the wave operator ap- 
pearing in equation (j8| but in the (id, fc)-domain. Its decomposition into real and imaginary 
parts highlights roles of propagation operator £ and dissipation operator T>, respectively, in 
that domain: 

Re|2(a;,fc) = (v% - v 2 s )k k T + (v 2 s k 2 - u 2 )lz (14) 
Im!D(u>,fc) = u(a 2 P — a%)k fc T + iDOsk 2 !^ (15) 

where Re and Im are the real and imaginary part operators. Then, the propagation is described 
by the real part of the dispersion matrix while the dissipation is described by its imaginary part. 

The relation of dispersion defined from the dispersion matrix by det SL(D,k) = 0, where 
det is the determinant operator, gives the propagation modes. More precisely, by introducing 
the dispersion manifold M as the set of all couples (w,fc) £ R x C 3 which satisfy the relation 
of dispersion, we show that (cj,jt) € M if and only if: 

2 2 

k 2 = k P (ui) := 5 or k 2 = fc|(w) := 2 , , (16) 
v 2 p (d) v|(w) 

The dispersion manifold naturally appears as the union of two manifolds i.e. M = Mp U Ms 
which corresponds to P- waves and S-waves with complex velocities vp and Vs, respectively. 

We consider now the low attenuation case. From the definition of vp and vj this assumption 
corresponds to Dap/v% << 1 and Da 2 s /v% « 1. Then, from (|16|l we get: 

k P (u) » — ; fc s (w) sb — (17) 
VP vs 

Those approximations of fcp(ui) and fcs(c^) will be useful to derive expression of the Green 
function and correlation and also Ward identities. 



4 Elastic Green function. 

In this section, we introduce and compute in some representation domains the Green func- 
tion of an elastic medium with viscous damping. We discuss the possibility to extract medium 
parameters from this field. The elastic Green function derived holds for unbounded medium, 
however, it is still usable for bounded media. 

4.1 Green function in the (oj, A;)-domain. 

For an unbounded, isotropic, homogeneous solid media, the Green function depends only 
on the difference between its spatial parameters i.e. we can do the following substitution: 
G(t,x,x') o G(t,x-x'). 
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As the Green function is the response to an isotropic spatio-temporal impulsion source i.e 
f := 5(t,x)l3, we obtain from equation (|12p : 

l(w,fc)g(w,fc) = I| (18) 

In that domain and according to equation QlgJI. the Green function appears to be the inverse 
of the dispersion matrix i.e. G(ui, k) = 3D (a;, k)~ ■ Using expression (|13[l . we show that: 

^ = i^^ kk + w^^- k - k ) (19) 

=: G P (uj,k)kk T + G s (uj,k) (h ~kk T ) (20) 

where k := fc/||fc|| and ||.|| is the classical Euclidean norm. The equation above highlights 
the decoupling between P-waves and S-waves in the (cj, fc)-domain. The P-waves displacement 
field, of amplitude determined by Gp(u>, fc), is on the axis directed by k whereas S-waves dis- 
placements, of amplitude determined by Gs(w,fe), occur in the plane orthogonal to k. 

It is important to note that Gp(uj,k) is the Green function of a propagation medium of 
acoustic waves with viscous damping. We recall that equation ((9]) describes the propagation in 
such media. 



4.2 Green function in the (u, x)-domain. 

We compute now the solid Green function in the (cj, x)-domain. In that domain, fields are 
capped by a V. 

Applying inverse Fourier transform to equation ([19]) with respect to k, we show in appendix 
EH that: 

G(u, x) = e . n n r; ^ x x + A ,. ,, a , . I 3 — a; a; + O tt-W [h -3xx (21) 
G P (iv,x)xx T + G s {uj,x) (lg-xx T ^ +G p ,s{uj,x) (lg-3xx T ^ (22) 

where x := x/llsll an d O is the classical Landau notation for dominated functions. The 
complete expression of G is given in appendix 17.11 G can be decomposed into the sum of a 
near-field term Gp,s = 0(||£|| _2 ) and two far-field terms Gp and Gs = 0(||a;||). A coupling 
term between the two waves dominates in the expression of the near-field contribution: 

G{to,x) « Gj., S (w,i) (l 3 - 3tI T ) (23) 

The far-field term highlights a decoupling between the two types of waves: 

G(uj, x) ~ Gp(oj, x)x i T + Gs(cj,3;) ^3 — x (24) 

We can note that we only retrieve the far-field contribution by taking the trace i.e. Tr(G) = 
Gp + 2Gg. This property is true in all representation domains. 
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When low attenuation case is considered, the Green function can then be approximated in 
the (w,a:)-domain using (|21l) and (JTTJ) by: 



LjJ_ 



G(oj,x) 



4ir\\x 



Jp - -t , e *s / _ _3 
bo — 4i i dj V= — 



+ 



+ 



Mlgllg + i 



la 



3i i 



47r||x|| 3 o; 2 



(25) 



Equation (|25[) is nothing but the classical solid Green function [T] in the (uj, a;)-domain for ideal 
case where no dissipation occurs. This approximation gives an easy interpretable expression 
of the Green function. Indeed, we observe pure phases in the far-field contribution for the 
two types of waves which provide information on ||x||,t)p and vs. This information can be 
completed with the near-field term. 



4.3 Approximated Green function in the (t, zj-domain. 

Computation of a general expression for G(t,x) from equation (|21|l when dissipation occurs, 
is difficult. In the low attenuation case, we obtain: 

G(t x ) „ S{t-tp) T 8(t-t S ) / _~~T\_ tn |lf , lsl W / _ T\ , v 

where ip := ||£||/wp, ts ■= \\x\\/vs and Il [tp>ts] (t) = 1 if i e [tp,is] and il [tpits] (t) = 
otherwise. (|26|l is the classical solid Green function [T] for ideal case in the (t, a:)-domain. In 
this simplified situation, interpreting the far-field term is easy as we can observe two times of 
arrival corresponding to each type of waves. The particle displacement of P-waves occurs along 
the propagation axis and particle displacement of S-waves remains in the orthogonal plane. 
For the near-field contribution, the coupling term controls the behavior of G, and provides also 
information on and vs- 



5 Elastic Green correlation and Ward identities. 

In this section, we compute the solid Green correlation and we highlight its role in passive 
identification of solid media. Results are also expressed for the acoustic case and for the far-field 
and low attenuation case. Exact and approximated Ward identities are derived and compared 
to existing ones for a different damping model. 

5.1 Computations in the (u, A;)-domain. 

As the Green function is space-shift invariant, the Green correlation also satisfies this prop- 
erty i.e. we can do the following substitution: C(t,x,x') 4-¥ C(t, x — x'). 

Using equation (0, we can compute the solid Green correlation in the (cj, fc)-domain. In- 
deed, applying the Fourier transform to equation ([7]) and noting that G~ = G where f 
transpose and conjugate, we obtain: 

k) = £0 , fc)£(oj , k) ] (27) 

The relation between the solid Green correlation and the solid Green function is purely alge- 
braic in this domain. 
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We can deduce a Ward identity in the (w, fc)-domain using equation (|27[) an d applying 
the lemma established in appendix 17.21 with A = = G(w,fc) _1 which satisfies lemma 

assumptions because of the presence of dissipation: 

uC(w,fe) = -D(fc) _1 Im &(<*;,&) (28) 
where D(fc) := 6U ImS(a;, fe) satisfies: 

D(fc)" 1 = a P 2 k' 2 kk T + as 2 r 2 fa - M) (29) 

D is the dispersion matrix of the dissipation operator T>. Ward identity (|28[) shows that we 
can retrieve the imaginary part of the solid Green function from the solid Green correlation. 
In fact, by invoking Kramers- Kroni'g theorem (|28[) based on the time-causality of the Green 
function, we can retrieve the complete Green function from the Green correlation. The ma- 
tricial proportionality term D(fc) -1 highlights the fundamental role of dissipation to establish 
Q28p. In fact, ([28} is true for every dissipation operator which is a partial differential operator. 
When the dissipation operator is considered as a constant operator as in [6], [5] and [7], we 
retrieve that the Green correlation is directly proportional to the imaginary part of the Green 
function. In the solid case with viscous damping, the proportionally term (|29[) is matricial and 
each of its components is proportional to k~ 2 . 



We can now compute the solid Green correlation in the (cj, fc)-domain. Using Ward identity 
(I28[l and expression (|19[1 . we get: 

= : C P (uj,k)kk T + C s {Lu,k) (la -fcfc T ) (31) 
where C P (oj,fc) = |Gp(c<;,fc)| 2 and C S (u,k) = \G s {to,k)\ 2 . 



Expression (]30|) proves that retrieving the Green correlation is sufficient to estimate medium 
parameters: vp, vs, ap and as- More precisely, those parameters are contained in kp(to) and 
fc|(cj) which are the poles of C(ui, k). 



For P-waves, we obtain from (|28|) a Ward identity applicable to acoustic waves with viscous 
damping: 

wC P (u,fe) = -^Iin6p(wJ (32) 

In that domain, the acoustic Green correlation is proportional to the imaginary part of the 
acoustic Green function. According to equation, (f32|) dissipation appears to be fundamental 
even in this scalar case. The proportionally term l/(a%k 2 ) is inverse of the dispersion matrix 
(of size 1 x 1) of the operator a P A which is the dissipation operator appearing in equation ((9]). 

5.2 Computations in the (a;, x) -domain. 

We compute now the Green correlation in the (cj, i)-domain. Inverse Fourier transform of 
equation (|30[) with respect to k gives ( appendix 17.311 : 

— Cp(u, x)xx T + Cg(a;, x) (la. ~ xx T ^j + Cp.g(cj, x) — 3xx T ^j (33) 
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The complete expression of C is given in appendix 17.31 and shows that we can extract all the 
physical parameters from (|33p . The decomposition above is similar to the decomposition into 
near-field and far-field of G. As for the Green function, we can retrieve only the far-field con- 
tribution by taking the trace i.e. Tr(C) = Cp + 2Cs. This property is also true in all Fourier 
domains. 



With far-field and low attenuation assumptions, (|33[) and (|17|) leads to: 



V vp ) ~ x ~ x t | V vs L (i 3 x f (34) 

x\\u} 3 a% 47r||a;||a; 3 a| V= / 



C{uj,x) « , 



Relations (|33f) and (|34[) are fundamental as they prove that we can extract medium parameters 
from the Green correlation in the (a;, £)-domain. This domain is accessible from the recorded 
field by an array of sensors using classical methods of discrete Fourier transform as the FFT 
(Fast Fourier Transform) algorithm. Green correlation expression (|34[) shows that in far-field 
and low attenuation case, we can retrieve times of arrival tp — \\x\\/vp and ts = \ \x\\/ v s using 
a pseudo-pulsation estimator. 

From Ward identity (|28|) . we can derived an exact Ward identity in the (a;, x)-domain: 

ujg(uj,x) = -23 _1 Img(w,x) (35) 

where T>~ is the inverse of dissipation operator. This identity extends the classical one for 
acoustic waves in a sense that the Green correlation is "proportional" to the imaginary part 
of the Green function where the inverse of the dissipation appears in the proportionality term. 
However, in this solid case with viscous damping, the proportionality term is V^ 1 and it is not 
a constant operator. Then, (|35|l is not sufficient to retrieve directly the Green function from 
the Green correlation. 

Similarly, we show that in the acoustic case, the Ward identity is: 

ljCp(lo,x) = -ap 2 A _1 ImGp(tij,i) (36) 

The "proportionality" term is the inverse of the Laplacian operator which is the dissipation 
operator. Even in this scalar case, Ward identity (|36|) is still hard to interpret. 



Using relations (|25[) and (|34[) , we obtain the following approximated Ward identity: 

oj 3 Q(uj,x) w -^(x)- 1 Img(w l£ ) (37) 

where: 

dix)- 1 = ^Lxx T + ^ (h-xx T ) (38) 

Identity (|37[1 highlights a matricial proportionality between the imaginary part of the approxi- 
mated solid Green function and the approximated solid Green correlation. The proportionality 
term is co^ 3 d(x)^ 1 which form is interpretable using the dissipation term dT>/dt in propagation 
equation ((Sj and relations of dispersion (|17[) . Indeed, a;" 1 is the contribution of operator d/dt 
in the (lu, i)-domain. i^ _2 d(2;) _1 is the matricial contribution of T>, it has indeed a similar form 
as D(fc) -1 (|29p in the (w, x)-domain. The proportionality term k~ 2 in (|28[1 has become the 
proportionality term oj -2 in (|37[1 due to the approximated relations of dispersion (|17p which 
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link proportionality k 2 and uj 2 for the two modes. 



Conserving only P- waves contribution in equation (J37J) , it comes: 

2 

u} 3 Cp(uj,x) w -^f Im G P (u,i) (39) 



In that domain, this approximated Ward identity highlights a proportionality between acoustic 
Green function and correlation. 

5.3 Computations in the (t, a:)-domain. 

The solid Green correlation in the (t, i)-domain can be obtained by Fourier transforming 
equation (|33[) . The complicated analytical form of the exponential terms makes this calculation 
difficult. However, we can deduce from Ward identity (f35)l . an exact Ward identity in the (t, x)- 
domain: 

dC(t,x) 

m = -T> Odd £(*,£) (40) 

where OddG := 1/2 ^G — Q.J is the odd part of the Green function. In the acoustic case, 
Ward identity (|36[) leads to: 

9Cp £>& = -a^ a A- 1 OddGp(t,x) (41) 

As identities f|35p and (|36p . identities (|40p and (|4ip extend classical acoustic results [6], [5], [7] 
with constant dissipation, but are difficult to apply in practice because of the presence of the 
inverse of dissipation operator T> is the solid case and operator apA in the acoustic case. 

We can get over this last step in the far-field and low attenuation situation. From equation 
(|37|l . we obtain: 

d 3 C(t x) 

« -d(x) - 1 Odd G (t , x) (42) 

The presence of the third-time derivative is justified by the presence of the term uj 3 in identity 
(|37|l . Finally, in the solid case with viscous damping, when far-field propagation and low at- 
tenuation are considered, the third-time derivative of the solid Green correlation is matricially 
proportional to the odd part of the solid Green function. The proportionally term is due to the 
dissipation and (|37fl highlights its fundamental role to establish a Ward identity in that domain. 



Retaining only P-waves terms in approximated Ward identity (|42|) . we obtain: 

d' A C P (t,x) 



<9t3 



OddGp(J, £ ) (43) 



Again, we can easily make an analogy with the classical identities obtained in the acoustic 
framework when the dissipation is supposed to be constant [6], [5], [7]. In that case and 
domain, it is shown that the first-time derivative of the Green correlation is proportional to 
the odd part of the Green function with a proportionality term inversely proportional to the 
constant dissipation. For a viscous damping and for far-field and low attenuation framework, 
it is the third-time derivative which is directly proportional to the odd part of the Green 
function in the (t, a:)-domain. The difference derivative order comes only from the form of the 
dissipation. 
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6 Conclusion. 

Our contributions can be summarized into two points. Firstly, we introduced the Green 
correlation in a general context. We presented its fundamental role in passive identification in 
the sense that all medium parameters can be retrieved from this field and without retrieving 
the Green function from a Ward identity. 

Secondly, we applied this theory to acoustic and solid waves with viscous damping. We 
computed, in different representation domains, the Green function, the Green correlation and 
associated Ward identities. Those latter highlight the fundamental role played by the dissi- 
pation operator as its frequency law describes the proportionality between the elastic Green 
function and correlation. Equations (|33[1 an d (|34[) give Green correlation expressions in an 
accessible domain from field measurements. Those expressions provides all mechanical infor- 
mation on the medium which highlights Green correlation approach in passive identification of 
elastic media. We considered the far-field and low attenuation case, which provides interesting 
and interpretable Ward identities in (uj,x)- and (t, x)-domains. More precisely, in that case, it 
is the third time-derivative of the Green correlation which is proportional to the odd part of 
the Green function (|42[) . (|43[) . This is a direct consequence of viscous damping model and this 
result was compared to classical ones in acoustic with other damping models. 

The perspective pursued is to apply experimentally our results to solid media with an 
embedded 3-components instrumentation in order to retrieve the complete solid Green function 
and correlation i. e. the whole matrices. It is a very interesting perspective as it has never been 
realised, from our knowledge. Role of sensor orientations will certainly be emphasized. 



All non direct calculus done to derive the Green function, the Green correlation and Ward 
identity are presented in this appendix. 

7.1 Computations for the elastic Green function. 

First, we need the following lemma: 



7 Appendix. 



Lemma. For a 6 C \ R, let Y(a, .) : k E K 3 
transform. We can show that: 



^ k' 2 -a 2 an< ^ i a ' ■) ^ e i nverse Fourier 




Then a i— » Y(a, ) can be prolonged to whole of C 



The starting point is (|19p . which can be rewritten as: 




(44) 
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Using the lemma, we have: 

fi(w,x) = a," 2 V \7 T Y{kp(oj),x) + uj- 2 (k 2 s (cj)X 3 - \7\/ T ) Y {k s (uj),x) 

= 7 I, I, 5 — m + ; I, I, , — (13-XX + 

[(i||*||M") - 1) e-^N^) + (-iMksH + 1) e -^ k ^] = ^2 

This equation gives the complete and exact expression of elastic Green function with viscous 
damping in the (cj, x)-domain. 

7.2 Fundamental lemma for Ward identity. 

We consider the following lemma: 

Lemma. Let A be an invertible matrix such that ImA := 1/2(A — A*) is invertible and 
AA f = A 1 A. Then we have: A~ 1 A t ~ 1 = - flmAV 1 Im ( A' 1 



Proof. The lemma is a direct consequence of the equality: — Im ( A 

This equality can be proved by expanding the left hand side using the definition of Im A, 

7.3 Computations for the elastic Green correlation. 

Using equations 1)270 and ()44)l , we have: 



then, 
C(ui, x) 



[kp 2 {Lo)Y (k P (oj),k)kk T ] Im [Y {k s (uj),k) (h- k s 2 (uj)kk T ^ 



Im [k P 2 (uj) V V T Y (k P (oj),x)] Im [(l 3 - fc- 2 (u,) V V T ) Y (k s (uj),x)] 



' A rvp u & Ct 2 g 



Im (V'H^H-A Im ( e -i\\*\\ks(u>)) 
V J_~~t V J_ ( l3 

47r||a;||a; 3 ap 4m\ \x\ \uj 3 a% \= 



47t||k| \bJ 3 Ctg 

rillTllfcr.^.! 

Im 



jjjg[|fcp(g) - l e -i\\x\\k P (w) , -i||£||fcs(aj) + 1 c -i\\x\\k s (u) 
a 2 p k P (to) a|fe|(a;) 



I3 — 3x x 

47r||x|| 3 tJ 3 



This is the complete and exact expression of the elastic Green correlation with viscous damping 
in the (ui, x)-domain. 
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Abstract 

This paper presents a new approach on passive identification of elastic propagation media. 
As passive identification relies upon noise field correlation, an original and perhaps more natural 
approach is developed which consists in considering the Green correlation. It is shown that 
Green correlation contains all medium parameters and provides an appealing alternative to 
the classical Green function estimation through a Ward identity. Our formalism allows to 
extend classical scalar passive identification models to vectorial ones and to take into account 
realistic dissipation models. This approach is applied to acoustic and solid waves with viscous 
damping. 

1 Introduction 

Passive identification of a propagation medium consists in retrieving medium parameters, 
by using uncontrolled noise fluctuations [14]. Such an idea has long been pursued in acoustics 
([10], [14]) and seismology ([3], [15]) and gave rise to numerous applications and experimental 
validations ([13], [TTJ, [6]). 

Preceding studies ([10], [3], [13], [TTJ, [6]) rely upon the estimation of the Green function of 
the medium. Such estimation is made possible by exploiting a Ward identity [16], which relates 
the noise correlation function to the Green function ([10], [14], [5]). The fundamental role of 
dissipation in Ward identity was outlined in [6], were dissipation is assumed to be constant; 
however, a constant dissipation model is hardly acceptable from a physical point of view ([8], 
|12j). and needs to be further discussed. 

The contribution of the paper is twofold. Firstly, an alternative approach to passive iden- 
tification is proposed, which is based upon the Green correlation function. Green correlation 
is introduced as the correlation of a propagated white noise field [7]. We show that Green 
correlation contains all physical parameters that need to be identified for a complete character- 
ization of the propagation medium. The motivation for introducing Green correlation comes 
from classical system identification theory, where noise based identification relies strongly on 
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the transformation of second order statistics through linear systems [9] . 

Secondly, we highlight the role of Green correlation in the framework of elastic waves with 
viscous damping. Green function and Green correlation are explicitly formulated in the time 
and space Fourier formalism. Ward identities are also developed. Furthermore, for low viscous 
damping and far-field assumptions, a new formulation of the Ward identity involving third 
order time derivative, is derived. 

The organization of the paper is the following. In section 2, we introduce passive iden- 
tification through a linear system approach. We recall the definition of the Green function 
and its role in medium identification. Cross-correlation of random fields is recalled, and white 
noise notion is introduced. Then, we define the Green correlation, we show its role in passive 
identification and we discuss its relation with the Green function through Ward identities. 

In section 3, solid waves equation with viscous damping is presented. Dispersion matrix 
and relation of dispersion are introduced. Acoustic waves propagation is seen as a particular 
case of solid propagation study. 

Solid and acoustic Green functions are computed in section 4. We highlight their role 
in solid and acoustic medium identification, respectively. Approximation of far-field and low 
attenuation is also considered. This particular, but realistic, case allows to derive explicit ex- 
pressions. 

Solid and acoustic Green correlation are computed in section 5. The role of Green correla- 
tion in passive identification is emphasized in those practical cases. Ward identities are derived 
and compared to existing ones for a constant damping model. 

2 Green function and Green correlation of a linear 
propagation medium 

In this section, we adopt a linear system approach to describe medium identification. We 
introduce the Green function and we relate it to the classical impulsional response of linear 
systems. White noise is introduced and its existence is discussed. From white noise, we define 
the Green correlation and we highlight its natural role in passive identification. 

2.1 Medium and fields: a system approach. 

We denote by u(t, x) the value of field u at time t and position x. When u has one compo- 
nent it is said scalar, otherwise it is said vectorial. 

A propagation medium can be seen as a system where the source field f is the input and 
the generated field u is the output of the system. Those two fields are related by a relation of 
the following type: 

u=<9(f) (1) 

where operator & is a representation of medium properties: physical laws, boundary and initial 
conditions. The medium is said linear when is a linear operator. Only linear medium are 
considered in this paper. 

With that formalism, medium identification can be seen as system identification where 
parameters are physical, like attenuation or propagation speed, and geometrical, like distance 
or orientation between sensors. 
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2.2 Green function of a linear medium. 



A linear medium X satisfies the superposition theorem i.e. the value u(t,x) of the generated 
field can be seen as the superposition of all contributions of elementary sources f [if , x[)dt' dx[ 
emitted at time t' during dt' period in the volume centered in x' and of dimensions da:', for all 
times t' and points x'. Mathematically, this can be written u(t, x) — f R .. x G(t, x, t\ x')f(t', x')dt'dx' 
and simplified by introducing the generalised convolution £g>T,s as: 



G is called the Green function of the medium as it is the kernel of the operator appearing 
in equation ([TJ. All medium mechanical parameters are contained in its expression. This high- 
lights the importance of retrieving the Green function in medium identification. 

Physically, the i-th column of G corresponds to the medium response of a spatio-temporal 
impulsion directed by the i-th axis of the reference. According to this interpretation, G is 
sometimes called " impulsional response" of the medium in reference to the classical impulsional 
response of a linear system [9]. Then, active identification consists in emitting spatio-temporal 
impulsions to retrieve the Green function and then to estimate model parameters [9]. 

We consider only time-shift invariant media. This property implies that G depends only 
on the times difference appearing in its parameters i.e. we can do the following substitution: 
G{t,x,t',x') o G(t - t',x,x'). 

2.3 Cross-correlation of stochastic fields. 

In passive identification, source fields are not controlled. The principle relies on record- 
ing noise sources and using their statistical properties to retrieve medium parameters. With 
stochastic source fields, the analysis has to be performed from the cross-correlation of the 
generated field u defined as: 



where E and T are the expectation operator and the transposition operator, respectively. 

We consider only stationary fields, this assumption is not a strong constraint in practice. 
In that case, cross-correlation depends only on the times difference appearing in its parameters 
i.e. we can do the following substitution: C(t,x,t' ,x') <-> C(t — t',x,x'). Furthermore, when 
fields are also considered ergodic in time, cross-correlation can be computed using the following 
formula: 



This formula is fundamental to approximate the cross-correlation of a field recorded by an 
array of sensors. 

2.4 White noise. 

By definition, a white noise is a field which value of a given time, position and direction is 
uncorrelated to any other value taken at all other times, positions and directions. Mathemati- 
cally, the cross-correlation of a such field f is a spatio-temporal isotropic impulsion: 



u = G <g>T,S f 



(2) 




(3) 




(4) 



C { (t,x, t' , x) — 5(t,t')5(x,x')I 



(5) 
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where p is the number of components of f , I is the p x p identity matrix and 5 is the Dirac 
distribution. 

A white noise has no physical reality in a sense that it has an infinite power. However, 
in practice the temporal whiteness is only needed in a limited frequency band. This latter is 
defined by the used instrumentation. The classical approach [3] to justify that ambient noise 
converges to a source with a spatio-isotropic whiteness consists to see the medium as a chaotic 
dynamical system. Then, according to equipartition theorem, it exists a time after which a 
coherent source snared in the medium becomes spatially and isotropically white. This time, 
called mixing time, depends on the frequency band, medium geometry and heterogeneity. With 
those considerations, cross-correlation of ambient noise is stacked during a sufficient long time 
in order to obtain a contribution of an approximated spatio-isotropic noise [3J, [BJ. 

2.5 Green correlation. 

As the Green function is the field generated by a spatio-temporal isotropic source i.e. 
f (t, x) = 5(t, t')8(x, , we define by analogy the Green correlation C by the cross-correlation 
of a field generated by a white noise source. This function, introduced in [7], plays by definition 
a fundamental role in passive identification. 

We can precise Green correlation expression using definition ((3J and equation ([2]). Indeed, 
for every generated field u, C can be expressed as: 

C u =g® T ,sC ( ® T| sG" (6) 

where G~ (t,x,x') := G(— t, x, x') T . It is important to note that to establish j6]), we use the 
property: G(t,x,x') = G(t, x , x), true for all times t and all couples of positions (x,x'), 
according to spatial reciprocity. Equation (J6j) is the "order two" version of equation ((2}. 
When the source is a white noise, we obtain by using equation |[6j) an expression of the Green 
correlation: 

C:=G® T ,sGr (7) 

This shows the fundamental importance of the Green correlation in medium identification when 
statistical properties of ambient sources are taking into account. Equations ([6j and (J7J show 
that "perfect" white noise is for passive identification what "perfect" impulsion is for active 
identification. 

Green correlation does not appear in the literature as a fundamental field to retrieve in 
order to estimate medium parameters. Generally, the cross-correlation of a field generated by 
a white noise through a propagation medium is directly related to the Green function by the 
authors. This relation is called Ward identity and relates actually the Green correlation to the 
Green function. 

3 Elastic propagation with viscous damping. 

In this section, we recall acoustic and solid propagation with viscous damping equations 
|12| . Classical results are derived with a vectorial formalism which is useful to describe solid 
waves propagation, in particular the coupling between the P-waves and the S- waves. Dispersion 
matrix and relation of dispersion are presented in order to compute the acoustic and solid Green 
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function, the acoustic and solid Green correlation and Ward identities. Low attenuation case 
is also studied because it provides physical interpretations. 



3.1 Solid waves equation with viscous damping. 

We consider an elastic, homogeneous, isotropic and linear solid medium. Let (A, fi) be the 
Lame parameters expressed in N.m~ 2 , (x,v) be the viscous damping parameters for P- waves 
and S- waves, respectively, expressed in N.s.m~ 2 , and p be the density of the medium expressed 
in kg.m" 3 . Let f be a 3-components causal spatio-temporal displacement source and u the 
3-components displacement field. From Newton theorem and Hooke law, we obtain the solid 
waves equation with viscous damping [12] : 



I3 is the identity operator of 3-components fields, 



u = f (8) 



C := — - — VV T — — AI3 (Propagation operator) 
P P 

V := - X + r, X7X7 T - -AJ 3 (Dissipation operator) 
P P 

where A is the Laplacian operator and V is the gradient operator. 



According to the Helmholtz- Hodge theorem [5] , each displacement field u satisfying equation 
{HJ when there is no source i.e. f = 0, can be decomposed asu = V(/>p + V A ^> , respectively, 
where A is the cross product. <f>p is a scalar field, and, ip is a two-components field satisfying: 

4>p = (9) 

t s =° ( 10 ) 

with vp := ^j^, v% := E , a.% := 2t ^ 2 and a| := 2 . This shows that a solid wave is the 
contribution of two modes satisfying the classical acoustic waves equation with viscous damping 
|12j . <j)p is a pressure wave, called P-wave, of speed vp and attenuation ap. ijj is a shear 
waves, called 5- wave, of speed vs and attenuation as- 



9 2 9 A 2 A 

__ Qp _ A _ VpA 



° T 2 A 2 A 



3.2 Dispersion matrix and relation of dispersion. 

For a deterministic field u, we denoted by: 

u(w,fc) := / n(^)-£) e_1 ^ t_ ~ ~^dtdx (11) 

Jmxm 3 

its Fourier transform in the (w, fc)-domain where ui and k are the frequency variables associated 
with t and x, respectively. 

Applying time and space Fourier transform to equation (|8]), we obtain the algebraic relation: 

S(w,fc)u(w,fe) = ffw.fc) (12) 
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where is the dispersion matrix: 

g(w,fc) := (v%(u) - v|(cj))fcfc T + (v|(o;)A: 2 - w 2 )I| (13) 

with: Vp(oj) := up + iujap, v|(w) := u| +iwa|, k 2 := fc T fc and I3 is the 3x3 identity matrix. 

The dispersion matrix is fundamental as it plays the same role as the wave operator ap- 
pearing in equation (j8| but in the (id, fc)-domain. Its decomposition into real and imaginary 
parts highlights roles of propagation operator £ and dissipation operator T>, respectively, in 
that domain: 

Re|2(a;,fc) = (v% - v 2 s )k k T + (v 2 s k 2 - u 2 )lz (14) 
Im!D(u>,fc) = u(a 2 P — a%)k fc T + iDOsk 2 !^ (15) 

where Re and Im are the real and imaginary part operators. Then, the propagation is described 
by the real part of the dispersion matrix while the dissipation is described by its imaginary part. 

The relation of dispersion defined from the dispersion matrix by det SL(D,k) = 0, where 
det is the determinant operator, gives the propagation modes. More precisely, by introducing 
the dispersion manifold M as the set of all couples (w,fc) £ R x C 3 which satisfy the relation 
of dispersion, we show that (cj,jt) € M if and only if: 

2 2 

k 2 = k P (ui) := 5 or k 2 = fc|(w) := 2 , , (16) 
v 2 p (d) v|(w) 

The dispersion manifold naturally appears as the union of two manifolds i.e. M = Mp U Ms 
which corresponds to P- waves and S-waves with complex velocities vp and Vs, respectively. 

We consider now the low attenuation case. From the definition of vp and vj this assumption 
corresponds to Dap/v% << 1 and Da 2 s /v% « 1. Then, from (|16|l we get: 

k P (u) » — ; fc s (w) sb — (17) 
VP vs 

Those approximations of fcp(ui) and fcs(c^) will be useful to derive expression of the Green 
function and correlation and also Ward identities. 



4 Elastic Green function. 

In this section, we introduce and compute in some representation domains the Green func- 
tion of an elastic medium with viscous damping. We discuss the possibility to extract medium 
parameters from this field. The elastic Green function derived holds for unbounded medium; 
however, it is still usable for bounded media, as long as different bouncing waves may be inter- 
preted as coming from new sources. This will be valid if the boundaries are far enough from 
the points of interest. The same line of reasoning hold for dealing with potentially present 
weak heterogeneities. 

4.1 Green function in the (u, A;)-domain. 

For an unbounded, isotropic, homogeneous solid media, the Green function depends only 
on the difference between its spatial parameters i.e. we can do the following substitution: 
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G(t,x,x') •*->• G(t,x-x'). 



As the Green function is the response to an isotropic spatio-temporal impulsion source i.e 
f := 8(t, x)l3_, we obtain from equation (|12[) : 

S(w,fe)G(w,fe) = Is (18) 



In that domain and according to equation (|18|l . the Green function appears to be the inverse 
of the dispersion matrix i.e. G(w, fc) = ®.{lo, fc) . Using expression (|13|l . we show that: 

=: Gp(Lj,fc)fcfc T + G s (a;,fc) (is -fcfc T ) (20) 



where fc := fc/||fc|| and ||.|| is the classical Euclidean norm. The equation above highlights 
the decoupling between P-waves and S-waves in the (cj, fc)-domain. The P-waves displacement 
field, of amplitude determined by Gp(uj,k), is on the axis directed by k whereas S-waves dis- 
placements, of amplitude determined by Gs(oj, k), occur in the plane orthogonal to k. 

It is important to note that Gp(cj,k) is the Green function of a propagation medium of 
acoustic waves with viscous damping. We recall that equation ((9]) describes the propagation in 
such media. 

4.2 Green function in the (u, x)-domain. 

We compute now the solid Green function in the (lj, s)-domain. In that domain, fields are 
capped by a V. 



Applying inverse Fourier transform to equation (|19l) with respect to k, we show in appendix 
imthat: 

&m= e ^^^ T + e ^m^ (2D 



47T I \ x 1 1 Vp (oj) 4ty \ \ x\ I (uj) 

G P (oj,x)xx T + G s {uj,x) (l|-ii T ) +G p , s {uj,x) (lg-3xx T ^ (22) 

where x_ := and O is the classical Landau notation for dominated functions. The 

complete expression of G is given in appendix 17.11 G can be decomposed into the sum of a 
near-field term Gp,s = 0(||3l||~ 2 ) and two far-field terms Gp and Gs = 0(||a<ll)- A coupling 
term between the two waves dominates in the expression of the near-field contribution: 

Q{ui,x) « Gp, s (oj,x) (l 3 - 3xx T ") (23) 

The far-field term highlights a decoupling between the two types of waves: 

G(u, x) ~ Gp(uj,x)xx T + Gs(w,2)) (I3 — x i T ) (24) 



We can note that we only retrieve the far-field contribution by taking the trace i.e. Tr(G) = 
Gp + 2Gs- This property is true in all representation domains. 
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When low attenuation case is considered, the Green function can then be approximated in 
the (w,a:)-domain using (|21l) and (JTTJ) by: 



LjJ_ 



G(oj,x) 



4ir\\x 



Jp - -t , e *s / _ _3 
bo — 4i i dj V= — 



+ 



+ 



Mlgllg + i 



la 



3i i 



47r||x|| 3 o; 2 



(25) 



Equation (|25[) is nothing but the classical solid Green function [T] in the (uj, a;)-domain for ideal 
case where no dissipation occurs. This approximation gives an easy interpretable expression 
of the Green function. Indeed, we observe pure phases in the far-field contribution for the 
two types of waves which provide information on ||x||,t)p and vs. This information can be 
completed with the near-field term. 



4.3 Approximated Green function in the (t, zj-domain. 

Computation of a general expression for G(t,x) from equation (|21|l when dissipation occurs, 
is difficult. In the low attenuation case, we obtain: 

G(t x ) „ S{t-tp) T 8(t-t S ) / _~~T\_ tn |lf , lsl W / _ T\ , v 

where ip := ||£||/wp, ts ■= \\x\\/vs and Il [tp>ts] (t) = 1 if i e [tp,is] and il [tpits] (t) = 
otherwise. (|26|l is the classical solid Green function [T] for ideal case in the (t, a:)-domain. In 
this simplified situation, interpreting the far-field term is easy as we can observe two times of 
arrival corresponding to each type of waves. The particle displacement of P-waves occurs along 
the propagation axis and particle displacement of S-waves remains in the orthogonal plane. 
For the near-field contribution, the coupling term controls the behavior of G, and provides also 
information on and vs- 



5 Elastic Green correlation and Ward identities. 

In this section, we compute the solid Green correlation and we highlight its role in passive 
identification of solid media. Results are also expressed for the acoustic case and for the far-field 
and low attenuation case. Exact and approximated Ward identities are derived and compared 
to existing ones for a different damping model. 

5.1 Computations in the (u, A;)-domain. 

As the Green function is space-shift invariant, the Green correlation also satisfies this prop- 
erty i.e. we can do the following substitution: C(t,x,x') 4-¥ C(t, x — x'). 

Using equation (0, we can compute the solid Green correlation in the (cj, fc)-domain. In- 
deed, applying the Fourier transform to equation ([7]) and noting that G~ = G where f 
transpose and conjugate, we obtain: 

k) = £0 , fc)£(oj , k) ] (27) 

The relation between the solid Green correlation and the solid Green function is purely alge- 
braic in this domain. 



8 



We can deduce a Ward identity in the (w, fc)-domain using equation (|27[) an d applying 
the lemma established in appendix 17.21 with A = = G(w,fc) _1 which satisfies lemma 

assumptions because of the presence of dissipation: 

uC(w,fe) = -D(fc) _1 Im &(<*;,&) (28) 
where D(fc) := 6U ImS(a;, fe) satisfies: 

D(fc)" 1 = a P 2 k' 2 kk T + as 2 r 2 fa - M) (29) 

D is the dispersion matrix of the dissipation operator T>. Ward identity (|28[) shows that we 
can retrieve the imaginary part of the solid Green function from the solid Green correlation. 
In fact, by invoking Kramers- Kroni'g theorem (|28[) based on the time-causality of the Green 
function, we can retrieve the complete Green function from the Green correlation. The ma- 
tricial proportionality term D(fc) -1 highlights the fundamental role of dissipation to establish 
Q28p. In fact, ([28} is true for every dissipation operator which is a partial differential operator. 
When the dissipation operator is considered as a constant operator as in [6], [5] and [7], we 
retrieve that the Green correlation is directly proportional to the imaginary part of the Green 
function. In the solid case with viscous damping, the proportionally term (|29[) is matricial and 
each of its components is proportional to k~ 2 . 



We can now compute the solid Green correlation in the (cj, fc)-domain. Using Ward identity 
(I28[l and expression (|19[1 . we get: 

= : C P (uj,k)kk T + C s {Lu,k) (la -fcfc T ) (31) 
where C P (oj,fc) = |Gp(c<;,fc)| 2 and C S (u,k) = \G s {to,k)\ 2 . 



Expression (]30|) proves that retrieving the Green correlation is sufficient to estimate medium 
parameters: vp, vs, ap and as- More precisely, those parameters are contained in kp(to) and 
fc|(cj) which are the poles of C(ui, k). 



For P-waves, we obtain from (|28|) a Ward identity applicable to acoustic waves with viscous 
damping: 

wC P (u,fe) = -^Iin6p(wJ (32) 

In that domain, the acoustic Green correlation is proportional to the imaginary part of the 
acoustic Green function. According to equation, (f32|) dissipation appears to be fundamental 
even in this scalar case. The proportionally term l/(a%k 2 ) is inverse of the dispersion matrix 
(of size 1 x 1) of the operator a P A which is the dissipation operator appearing in equation ((9]). 

5.2 Computations in the (a;, x) -domain. 

We compute now the Green correlation in the (cj, i)-domain. Inverse Fourier transform of 
equation (|30[) with respect to k gives ( appendix 17.311 : 

— Cp(u, x)xx T + Cg(a;, x) (la. ~ xx T ^j + Cp.g(cj, x) — 3xx T ^j (33) 
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The complete expression of C is given in appendix 17.31 and shows that we can extract all the 
physical parameters from (|33p . The decomposition above is similar to the decomposition into 
near-field and far-field of G. As for the Green function, we can retrieve only the far-field con- 
tribution by taking the trace i.e. Tr(C) = Cp + 2Cs. This property is also true in all Fourier 
domains. 



With far-field and low attenuation assumptions, (|33[) and (|17|) leads to: 



V vp ) ~ x ~ x t | V vs L (i 3 x f (34) 

x\\u} 3 a% 47r||a;||a; 3 a| V= / 



C{uj,x) « , 



Relations (|33f) and (|34[) are fundamental as they prove that we can extract medium parameters 
from the Green correlation in the (a;, £)-domain. This domain is accessible from the recorded 
field by an array of sensors using classical methods of discrete Fourier transform as the FFT 
(Fast Fourier Transform) algorithm. Green correlation expression (|34[) shows that in far-field 
and low attenuation case, we can retrieve times of arrival tp — \\x\\/vp and ts = \ \x\\/ v s using 
a pseudo-pulsation estimator. 

From Ward identity (|28|) . we can derived an exact Ward identity in the (a;, x)-domain: 

ujg(uj,x) = -23 _1 Img(w,x) (35) 

where T>~ is the inverse of dissipation operator. This identity extends the classical one for 
acoustic waves in a sense that the Green correlation is "proportional" to the imaginary part 
of the Green function where the inverse of the dissipation appears in the proportionality term. 
However, in this solid case with viscous damping, the proportionality term is V^ 1 and it is not 
a constant operator. Then, (|35|l is not sufficient to retrieve directly the Green function from 
the Green correlation. 

Similarly, we show that in the acoustic case, the Ward identity is: 

ljCp(lo,x) = -ap 2 A _1 ImGp(tij,i) (36) 

The "proportionality" term is the inverse of the Laplacian operator which is the dissipation 
operator. Even in this scalar case, Ward identity (|36|) is still hard to interpret. 



Using relations (|25[) and (|34[) , we obtain the following approximated Ward identity: 

oj 3 Q(uj,x) w -^(x)- 1 Img(w l£ ) (37) 

where: 

dix)- 1 = ^Lxx T + ^ (h-xx T ) (38) 

Identity (|37[1 highlights a matricial proportionality between the imaginary part of the approxi- 
mated solid Green function and the approximated solid Green correlation. The proportionality 
term is co^ 3 d(x)^ 1 which form is interpretable using the dissipation term dT>/dt in propagation 
equation ((Sj and relations of dispersion (|17[) . Indeed, a;" 1 is the contribution of operator d/dt 
in the (lu, i)-domain. i^ _2 d(2;) _1 is the matricial contribution of T>, it has indeed a similar form 
as D(fc) -1 (|29p in the (w, x)-domain. The proportionality term k~ 2 in (|28[1 has become the 
proportionality term oj -2 in (|37[1 due to the approximated relations of dispersion (|17p which 
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link proportionality k 2 and uj 2 for the two modes. 



Conserving only P- waves contribution in equation (J37J) , it comes: 

2 

u} 3 Cp(uj,x) w -^f Im G P (u,i) (39) 



In that domain, this approximated Ward identity highlights a proportionality between acoustic 
Green function and correlation. 

5.3 Computations in the (t, a:)-domain. 

The solid Green correlation in the (t, i)-domain can be obtained by Fourier transforming 
equation (|33[) . The complicated analytical form of the exponential terms makes this calculation 
difficult. However, we can deduce from Ward identity (f35)l . an exact Ward identity in the (t, x)- 
domain: 

dC(t,x) 

m = -T> Odd £(*,£) (40) 

where OddG := 1/2 ^G — Q.J is the odd part of the Green function. In the acoustic case, 
Ward identity (|36[) leads to: 

9Cp £>& = -a^ a A- 1 OddGp(t,x) (41) 

As identities f|35p and (|36p . identities (|40p and (|4ip extend classical acoustic results [6], [5], [7] 
with constant dissipation, but are difficult to apply in practice because of the presence of the 
inverse of dissipation operator T> is the solid case and operator apA in the acoustic case. 

We can get over this last step in the far-field and low attenuation situation. From equation 
(|37|l . we obtain: 

d 3 C(t x) 

« -d(x) - 1 Odd G (t , x) (42) 

The presence of the third-time derivative is justified by the presence of the term uj 3 in identity 
(|37|l . Finally, in the solid case with viscous damping, when far-field propagation and low at- 
tenuation are considered, the third-time derivative of the solid Green correlation is matricially 
proportional to the odd part of the solid Green function. The proportionally term is due to the 
dissipation and (|37fl highlights its fundamental role to establish a Ward identity in that domain. 



Retaining only P-waves terms in approximated Ward identity (|42|) . we obtain: 

d' A C P (t,x) 



<9t3 



OddGp(J, £ ) (43) 



Again, we can easily make an analogy with the classical identities obtained in the acoustic 
framework when the dissipation is supposed to be constant [6], [5], [7]. In that case and 
domain, it is shown that the first-time derivative of the Green correlation is proportional to 
the odd part of the Green function with a proportionality term inversely proportional to the 
constant dissipation. For a viscous damping and for far-field and low attenuation framework, 
it is the third-time derivative which is directly proportional to the odd part of the Green 
function in the (t, a:)-domain. The difference derivative order comes only from the form of the 
dissipation. 
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6 Conclusion. 

Our contributions can be summarized into two points. Firstly, we introduced the Green 
correlation in a general context. We presented its fundamental role in passive identification in 
the sense that all medium parameters can be retrieved from this field and without retrieving 
the Green function from a Ward identity. 

Secondly, we applied this theory to acoustic and solid waves with viscous damping. We 
computed, in different representation domains, the Green function, the Green correlation and 
associated Ward identities. Those latter highlight the fundamental role played by the dissi- 
pation operator as its frequency law describes the proportionality between the elastic Green 
function and correlation. Equations (|33[1 an d (|34[) give Green correlation expressions in an 
accessible domain from field measurements. Those expressions provides all mechanical infor- 
mation on the medium which highlights Green correlation approach in passive identification of 
elastic media. We considered the far-field and low attenuation case, which provides interesting 
and interpretable Ward identities in (uj,x)- and (t, x)-domains. More precisely, in that case, it 
is the third time-derivative of the Green correlation which is proportional to the odd part of 
the Green function (|42[) . (|43[) . This is a direct consequence of viscous damping model and this 
result was compared to classical ones in acoustic with other damping models. 

The perspective pursued is to apply experimentally our results to solid media with an 
embedded 3-components instrumentation in order to retrieve the complete solid Green function 
and correlation i. e. the whole matrices. It is a very interesting perspective as it has never been 
realised, from our knowledge. Role of sensor orientations will certainly be emphasized. 



All non direct calculus done to derive the Green function, the Green correlation and Ward 
identity are presented in this appendix. 

7.1 Computations for the elastic Green function. 

First, we need the following lemma: 



7 Appendix. 



Lemma. For a 6 C \ R, let Y(a, .) : k E K 3 
transform. We can show that: 



^ k' 2 -a 2 an< ^ i a ' ■) ^ e i nverse Fourier 




Then a i— » Y(a, ) can be prolonged to whole of C 



The starting point is (|19p . which can be rewritten as: 




(44) 
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Using the lemma, we have: 

fi(w,x) = a," 2 V \7 T Y{kp(oj),x) + uj- 2 (k 2 s (cj)X 3 - \7\/ T ) Y {k s (uj),x) 

= 7 I, I, 5 — m + ; I, I, , — (13-XX + 

[(i||*||M") - 1) e-^N^) + (-iMksH + 1) e -^ k ^] = ^2 

This equation gives the complete and exact expression of elastic Green function with viscous 
damping in the (cj, x)-domain. 

7.2 Fundamental lemma for Ward identity. 

We consider the following lemma: 

Lemma. Let A be an invertible matrix such that ImA := 1/2(A — A*) is invertible and 
AA f = A 1 A. Then we have: A~ 1 A t ~ 1 = - flmAV 1 Im ( A' 1 



Proof. The lemma is a direct consequence of the equality: — Im ( A 

This equality can be proved by expanding the left hand side using the definition of Im A, 

7.3 Computations for the elastic Green correlation. 

Using equations 1)270 and ()44)l , we have: 



then, 
C(ui, x) 



[kp 2 {Lo)Y (k P (oj),k)kk T ] Im [Y {k s (uj),k) (h- k s 2 (uj)kk T ^ 



Im [k P 2 (uj) V V T Y (k P (oj),x)] Im [(l 3 - fc- 2 (u,) V V T ) Y (k s (uj),x)] 



' A rvp u & Ct 2 g 



Im (V'H^H-A Im ( e -i\\*\\ks(u>)) 
V J_~~t V J_ ( l3 

47r||a;||a; 3 ap 4m\ \x\ \uj 3 a% \= 



47t||k| \bJ 3 Ctg 

rillTllfcr.^.! 

Im 



jjjg[|fcp(g) - l e -i\\x\\k P (w) , -i||£||fcs(aj) + 1 c -i\\x\\k s (u) 
a 2 p k P (to) a|fe|(a;) 



I3 — 3x x 

47r||x|| 3 tJ 3 



This is the complete and exact expression of the elastic Green correlation with viscous damping 
in the (ui, x)-domain. 
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